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Abstract. The s-ordered phase-sum and phase-difference distributions are consid- 
ered for Bell-like superpositions of two-mode coherent states. The distributions are 
sensitive, respectively, to the sum and difference of the phases of the entangled coher- 
ent states. They show loss of information about the entangled state and may take on 
negative values for some orderings s. 



PACS numbers: 42.50 Dv, 03.65 Ud 
1. Introduction 

The quantum phase issue is still an open problem, since there are several nonequivalent 
proposals to extend the notion of a phase from classical to quantum physics: construction 
of a phase operator canonically conjugate to the photon number operator 
derivation of a phase distribution by integrating a phase space distribution over the 
"radial" variables [SI E] or operational definition by a certain experimental set-up [Sj. 
Here we choose the phase space approach in order to discuss the phase properties of 
entangled two-mode coherent states. The two-mode (or bipartite) states have acquired a 
great interest due to their potential applications in secure communications [B], quantum 
teleportation |Hj as well as in questioning the fundamentals of quantum theory, e.g. 
in testing Bell inequalities :7J. On the other hand, the s-ordered multimode quasi- 
probability distributions can be sampled point-by-point using unbalanced homodyning 
[Hj. The corresponding phase distributions considered here then follow as marginal 
distributions. 

This paper is organized as follows. In Section|2]we calculate for the quasi-Bell state 
considered the s-ordered characteristic function and quasi-probability distribution. The 
corresponding s-ordered phase-sum and phase-difference distributions are derived and 
discussed in Section 01 The concluding Section 0] contains a brief summary. 

§ To whom correspondence should be addressed (aalexandrescu@tehfi.pub.ro). Now with Faculty of 
Electronics and Telecommunications - Optoelectronics Research Center, " Politehnica" University of 
Bucharest, RO-06f07f Bucharest, Romania. 
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2. Characteristic functions and quasi-probability distributions 

We consider two bosonic modes with annihilation (creation) operators a, b (a', b'). 
One-mode coherent states \a) (\/3)) are defined as eigenstates of a (b) with complex 
eigenvalues a = \a\ exp(i<p a ) (/? = |/3| exp(i(pp)). Denoting by \a, (3) = \a) (g> |/3) the 
two-mode coherent states, we consider quasi-Bell states of the form 

\*)=Ar(fi\a,l3)+u\-a,-13)), (1) 

where /i and v are complex parameters restricted by |/i| 2 + \v\ 2 = 1 and M is the 
normalization constant 

M = {l + 2Re(/xz/*) exp[-2(|a| 2 + \f3\ 2 )}}~ h . (2) 

The superposition (JTJ) represents a general two-mode Schrodinger cat state. Particular 
states obtained for (p,, u) = (l,±l)/\/2 and (fJL,v) = (1, ±i)/\/2 are referred to as 
even/odd Schrodinger cat and Yurke-Stoler e 9j states, respectively By a suitable 
choice of the parameters a, (3, p and v the state may exhibit maximal bipartite 
entanglement JU] or two-mode squeezing 

Given the density operator p of a system one defines the complex- valued s-ordered 
characteristic function by [TJl HB] : 

X (£, t;; s) = exp g(|£| 2 + |r/| 2 )) Tr{ Pj D(£, r,)}, (3) 

where D(£,rj) = exp(£at — £*a) exp^Jr — r/*&) is the two-mode displacement operator 
with complex displacements £ and rj, and s is a complex parameter. The discrete 
values s = 1,0,-1 correspond to normal, symmetric and antinormal ordering of the 
boson operators, respectively. The characteristic function of the state ((TJ with density 
operator p = |^} (\&| is easily obtained as 

17; s) = AT 2 exp (-i^i (|£| 2 + M 2 )) { I/,) 2 exp(£a* - T« + - V*P) 

+ |z/| 2 exp(-£a* + fa - n[3* + n* (3) 

+ exp(-2|a| 2 -2|/3| 2 ) [//Vexp(£a<* + fa + + ?7*/3) 

+ pv* exp(-£a* - fa - VP* ~ V*P)}} • (4) 

The s-ordered quasi-probability distribution can be derived from the characteristic 
function by taking its Fourier transform 

W(j, M = ^ / d 2 £ J d 2 V 77; s) exp( 7 f - 7*0 exp(5 V * - 5* V ). (5) 

In order for this transform to exist one must have Re(s) < 1. For s — — 1 one obtains 
the (positive valued) Q(or Husimi)-function, for s = the Wigner function and for s = 1 
the (highly singular) Glauber-Sudarshan P-representation. Inserting the characteristic 
function (J3J) in relation (jHJ) one obtains the following expression for the s-ordered quasi- 
probability distribution of the state (JIJ: 

W(j, 6; s) = ^^exp (-^(M 2 + |/?| 2 )) exp ("^(N 2 + l*| 2 )) 
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It has the property [W(7, 8; s)]* = W(j,8;s*) and thus is real-valued only for real 
values of s. In this case considered henceforth the distribution W(j, 8; s) has two major 
contributions: the first one containing \fi\ 2 , \v\ 2 and exponentials of real numbers is 
always positive (Gaussian terms), whereas the second one containing /iV, \w* and 
exponentials of imaginary numbers may take on negative values as well (interference 
terms). By varying the value of the parameter s one can enhance or suppress the 
interference terms in comparison with the Gaussian ones. This results in always positive 
distributions W(j, 8; s) only for s < —1. If — 1 < s < 1, the distributions take on 
negative values as well. 

3. Phase-sum and phase-difference distributions 

In experiments one measures only relative phases, e.g. with respect to the phase of a 
local oscillator as in homodyne detection . Therefore, it makes sense to address 

the problem of phase-difference and, complementary, of phase-sum distributions in the 
phase operator approach [TH] as well as in the phase space approach [T7] . 

The main steps in deriving the phase-sum and phase-difference distributions from 
the two-mode quasi-probability distribution (JHJ) are: 

(i) The normalization condition of W^, 8; s) is written in polar coordinates (7 = 
N exp(iy? 7 ), 8 = \8\ exp(i(p s )): 

roo roo rlix rlix 

\ d| 7 | / d|tf| / & V J d^| 7 ii<W7,M = l- (7) 
J Jo Jo Jo 

(ii) The phase variables ip 1 and ips are replaced by the phase-sum and -difference: 

0+ = if S + <Ay> 4>~ = fS ~ W ( 8 ) 

In going from two one-mode phases defined on 27r-domains to their sum and 
difference a problem occurs because + and 0_ have 47r-ranges according to (JHJ. 
This is not compatible with measurements that cannot tell a phase value (f) from 
cf) ± 2tt. To cast the phase-sum and -difference into 27r-ranges, we follow reference 
[T5] and replace W(j, 8; s) by 

W(| 7 |, \8\,(f>+, 0-! s)=- [W(bi, \8\,<p y , W, s) + W{\^\, \S\,(p^ + tt, ip 5 + vr; s)} (9) 

with (p 1 and (p$ on the right-hand side substituted in terms of <j)± from relation (jHJ). 
The quasi-probability distribution W is then normalized according to 

/*oo rliv r2n 

d| 7 | / d\8\ d0 + / d0_| 7 |HW(l7l,H,0+,0-;s) = i- (10) 



JO JO 
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(iii) The phase-sum (phase-difference) distribution is obtained by integration over the 
"radial" variables |7|, \6\ and over (</>+): 

/•oo roo r-l'K 

V {± \^s)= d| 7 | / d|*|/ d0 T | 7 ||5| W(| 7 |,|5U + ,0-; S ). (11) 
Jo Jo Jo 

Applying the above procedure to the quasi-probability distribution (jHJ), we first expand 
the exponential terms of the distribution JOJ) into series of Bessel functions then we 
integrate over the appropriate phase variable (</>_ or <p + ) and finally over the "radial" 
variables ^H] • We write the final result for the s-ordered phase-sum and phase-difference 
distributions in the form of a genuine phase distribution: 



V^(<P ± ;s) = - 

Z7T 



1 + 2E 



^'(s) cosn ' ^ 1 — ^ 



(12) 



n=l 

with state-dependent phases <p' ± = tpp ± (p a and coefficients c^(s) given by: 

4 ±) (*)=^ , |{4 +) (S)^ ) (S) 

+ (±)» 2 Re(^*) e- 2 (l«l 2 +l^ 2 ) 4")(g) JH(® } . (13) 
Here A/" is the normalization constant (j2J and we have introduced the notations 

/W(x) = v^e^(j^(x)±J„ ± i(z) > ) (14) 
V 2 2 / 

" f ^ T j(2x)te^M(f + l,n + l,±2x), (15) 

where J^(x) is the modified Bessel function of order v and M(a,b,x) is Kummer's 
confluent hypergeometric function. The relation (JT5|) follows from (|14|) by relating /^(a;) 
to M(a, 6, x) and using then one of the recurrence relations ^H]- When calculating 
the phase-sum and -difference distributions (fTTj) with <f>± spanning a 27r-domain, part 
of the information about the state is lost. This loss of information is due to the 
imposed 27r-periodicity in (j>± which renders the transformation (jHJ) nonbijective [T8J. In 
our case the superposition (JTJ) depends on two complex numbers /1 and v subject to 
the condition |/i| 2 + |z/| 2 = 1, so that only three real independent parameters are left. 
Correspondingly the density operator and therefore the quasi-probability distribution 
(JBJ) depend on the three parameters | | 2 — |z/| 2 , Re(/iz/*) and Im(/iz/*). In comparison, 
the phase distributions (|12|) depend only on Re(/iz/*). The terms containing |/i| 2 — \v\ 2 
and Im(//i/*) were lost already at the level of the quasi-probability distribution W 
when evaluated according to relation (JSJ). An additional loss of information occurs 
after performing the integration (JTTj) and is related with the dependence on the state- 
related phases tp a and (fp: they enter the distribution V^ + ' (V^) only through the 
corresponding combination ipp + ip a (ipp — ip a ). The right-hand side of ((H)) is a Fourier 
series expansion with the sine terms missing. This makes the distributions symmetric 
about the origine <j)± = 4>'±. 

We note that less information is lost when considering the marginal one-mode phase 
distributions. Thus, integrating W(7, 5; s) in relation (JHJ) over the remaining variables, 
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we obtain the phase distribution for mode 1 in the form 
if 00 

V^(^; s) = — 1 + 2 X; cos(2A# 7 - </?«)) 

Z7r I fe=i 

+ c&UOO cos((2A; - 1)(^ 7 - ip a )) + ^(s) sin((2fc - 1)(^ 7 - p a ))] } , (16) 
where the coefficients are given by 



c. 



(i) 

2fc-l 



= ^ 2 4 { 2 lm( ^ } e ~ 2(H2+i/3|2) 4ii(S }• (17) 

The marginal phase distribution V^(fs] s ) of mode 2 has the same form ()16|) in terms 
of the phase variable ip$ — tpp and corresponding coefficients c£\s) and d k 2 (s). The 
latter are obtained formally by interchanging |a| «-> |/3| in ()17jl. These one- mode phase 
distributions depend on all but one of the relevant parameters defining the state (fTj). 
Thus, the distribution of mode 1 (2) does not depend on the phase (pp (ip a ) of the other 
mode. Due to the sine terms present in (fTE|) . the distributions are no longer symmetric 
about the phase origine. Note also that the sine terms with even indices are still missing. 

In Figure [T] we display phase-difference and in Figure El phase-sum distributions 
for even (upper subfigures (a), (b)) and odd (lower subfigures (c), (d)) Schrodinger cat 
states. Subfigures (a) and (c) on the left show the phase distributions (fTi^ over the plane 
(|a| 2 , 4>± — 4>'±) for s = and |a| = \(3\. The dependence on the ordering parameter s and 
the phases <fi± — <p' ± is plotted for three values of s (-1, and 0.4) and for \a\ = |/3| = 1 
in subfigures (b) and (d) on the right. The phase distributions are symmetric about 
4>± = <fr± = (pp ± <f a , with a peak or dip at the symmetry point. For small values of |ck| 
they approach the uniform distribution as long as the normalization constant M in (JT3j) 
stays finite at |a| = 0. However, if 1 + 2Re(/iz/*) = 0, as in the example of Figure 1(c), 
a non-uniform distribution is approached. Regarding the dependence on the ordering 
parameter, the phase distributions corresponding to s > —1 may take on negative values 
and may develop strong oscillations, especially for s > 0. This is due to the fact that 
the interference terms in (jUJ) become more important and may overwhelm the Gaussian 
terms. In general, the smaller the value of the ordering parameter is, the smoother is the 
corresponding phase distribution. Note that quasi-probability distributions W taking 
on also negative values may nevertheless yield positive- valued phase distributions V, as 
can be seen from Figures 1(b) and 2(d) for s = 0. This is not surprising, as the marginal 
distributions of a quasi-distribution may be true distributions. 

Given the quasi-probability distribution (fT2*|) one can calculate the expectation value 
of various phase-dependent functions. Here we give the central trigonometric moments 
(we drop the indices ± on 0, cp' and c n ): 

(cosn(0 -4>')) = c n , (sinn(0-0')) = 0, (18) 
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the variances squared 

(A cosn(0 - 0')) 2 = \ (l - 24 + c 2 n) , (A sinn(0 - 0')) 2 = ^(1 - c 2n ) (19) 
and the mean value and variance squared of the phase for the 27r-window [0 O — 7r, o +7r]: 

00 / \n 

(0) = 0o + 2 Y, — c n sin n(0 o - 0') , (20) 

n=l n 

2 00 / \n 

(A0) 2 = - - ((0) - O ) 2 + 4 £ H-Cn cosn(0 o - 0'). (21) 

d n=l W 

The latter expressions simplify considerably by choosing 0o = 0' for the center of the 
window, i.e. 0o = 0'_j_ = y^g ± (p a for the phase-sum/phase-difference distributions. It 
follows from (fTH|) that the coefficients c n (s) are the expectation values of cos n(0— 0') and 
therefore must satisfy the bound |c n (s) | < 1 for a genuine phase distribution. This bound 
may be violated, however, for quasi-distributions taking on negative values. This is the 
reason why such phase distributions, as in Figure 2(d), may develop a (negative) dip 
instead of a (positive) peak at = 0'. Similar considerations apply to the coefficients 
c^(s) and d$(s), i = 1,2, of the one-mode phase distributions (fi7)|) . which are the 
expectation values of the corresponding cosine and sine functions multiplying them. 





Figure 1. Phase-difference distributions. Left: From the Wigner function (s = 0) for 
(a) even and (c) odd Schrodinger cat states with |a| = \/3\. Right: Dependence on the 

ordering parameter, s — — 1 ( ), s = (- -) and s — 0.4 ( ) for (b) even and 

(d) odd Schrodinger cat states with \a\ = \(3\ = 1. 




4. Summary 

In this paper we have considered various marginal phase distributions obtained from the 
s-ordered two-mode quasi-probability distribution for the entangled two-mode coherent 
state in Proper non-negative distributions are obtained for s < — 1 and improper 
or quasi-distributions for — 1 < s < 1. The latter take on negative values as well and 
may yield unphysical expectation values for the (central) trigonometric moments. The 
Glauber-Sudarshan P-representation corresponding to s — 1 is not defined. 

As expected and contrary to the quasi-probability distribution W(7, <5; s) in (0), the 
marginal phase distributions are not sensitive to the full set of parameters characterizing 
the state (JTJ). Thus, the phase distributions for one mode depend only on its own phase 
and do not feel the phase of the other mode. Similarly, the phase-sum and phase- 
difference distributions involve, respectively, only the sum ipp+ipa and difference ipp — ipa 
of the two state-related phases. In addition, the two latter distributions do not depend 
on |/i| 2 — |z/| 2 and Im(/xz/*), but only on Re(/ii/*), as a consequence of casting the 4ir- 
periodic distribution W (considered as a function of <p + and 0_) into a 27r-periodic one 
W according to relation (JHJ). Regarded as Fourier series expansions, the phase-sum and 
phase-difference distributions have only (even) cosine terms and are therefore symmetric 
about the origine. The one-mode phase distributions, however, are unsymmetric due to 
the (odd) sine terms. 
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